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Abstract. We perform a systematic theoretical and numerical study of 
anti-parallel two-dimensional magnetic reconnection with asymmetries in the 
plasma density and reconnecting magnetic held strength in addition to a bulk 
how shear across the reconnection site in the plane of the reconnecting helds, 
which commonly occurs at planetary magnetospheres. We analytically pre¬ 
dict the speed at which an isolated X-line is convected by the how, the re¬ 
connection rate, and the critical how speed at which reconnection no longer 
takes place for arbitrary reconnecting magnetic held strengths, densities, and 
upstream how speeds, and we conhrm the results with two-huid numerical 
simulations. The predictions and simulation results counter the prevailing 
model of reconnection at Earth’s dayside magnetopause which says recon¬ 
nection occurs with a stationary X-line for sub-Alfvenic magnetosheath how, 
reconnection occurs but the X-line convects for magnetosheath hows between 
the Alfven speed and double the Alfven speed, and reconnection does not 
occur for magnetosheath hows greater than double the Alfven speed. In par¬ 
ticular, we hnd that X-line motion is governed by momentum conservation 
from the upstream hows, which are weighted diherently in asymmetric sys¬ 
tems, so the X-line convects for generic conditions including sub-Alfvenic up¬ 
stream speeds. For the reconnection rate, as with symmetric reconnection, 
it drops with increasing how shear and there is a cutoh speed above which 
reconnection is not predominant. However, while the cutoh condition for sym¬ 
metric reconnection is that the diherence in hows on the two sides of the re¬ 
connection site is twice the Alfven speed, we hnd asymmetries cause the cut- 
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off speed for asymmetric reconnection to be higher than twice the asymmet¬ 
ric form of the Alfven speed. The stronger the asymmetries, the more the 
cutoff exceeds double the asymmetric Alfven speed. This is due to the fact 
that in asymmetric reconnection, the plasma with the smaller mass flux into 
the dissipation region contributes a smaller mass to the dissipation region, 
so the effect of its flow on opposing the release of energy by the reconnected 
magnetic helds is diminished and the reconnection is not suppressed to the 
extent previously thought. The results compare favorably with an observa¬ 
tion of reconnection at Earth’s polar cusps during a period of northward in¬ 
terplanetary magnetic held, where reconnection occurs despite the magne¬ 
tosheath how speed being more than twice the magnetosheath Alfven speed, 
the previously proposed suppression condition. These results are expected 
to be of broad importance for magnetospheric physics of Earth and other plan¬ 
ets; particular applications are discussed. 
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1. Introduction 

A key element controlling the interaction of the solar wind with Earth’s magnetosphere 
is the nature and efficiency of the magnetic reconnection process at the dayside mag¬ 
netopause. During reconnection, solar wind magnetic field lines effectively break and 
cross-connect with terrestrial magnetic field lines. As a consequence, solar wind plasma is 
able to enter the magnetosphere and the reconnected magnetic field lines convect tailward 
[Dungey, 1961]. This is a crucial aspect of space weather phenomena. The focus of this 
study is how the bulk flow of the solar wind around the magnetosphere itself affects the 
dayside reconnection process; we focus on (1) the convection of the reconnection site by 
the bulk flow, (2) the effect on the reconnection rate, and (3) the critical bulk flow speed 
above which reconnection is not the dominant effect. 

Bulk flow is expected to be most important when reconnection occurs near the polar 
cusps. This is because the bulk magnetosheath flow around the magnetopause acquires 
a potentially sizable component parallel or anti-parallel to the magnetic field when the 
reconnection site is near the cusps. This affects the reconnection site very differently than 
when reconnection is near the subsolar point, where the flow is predominantly not aligned 
with the reconnecting magnetosheath magnetic held. As pointed out by Dungey [1963], 
magnetic reconnection is likely to occur near the polar cusps when the interplanetary mag¬ 
netic held (IMF) has a northward component. Detections of high latitude reconnection are 
quite common [Gosling et ai, 1986, 1991, 1996; Kessel et ai, 1996; Fuselier et ai, 2000a; 
Onsager et ai, 2001; Avanov et ai, 2001; Federov et ai, 2001; Phan et ai, 2001; Frey 
et al, 2003; Phan et a/., 2003; Lavraud et al, 2004, 2005; Retino et ai, 2005; Retino et al, 
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2006; Phan et al, 2006, 2007; Hasegawa et al, 2008; Fuselier et ai, 2010, 2012, 2014a; 
Muzamil et al, 2014; Wilder et al., 2014], 

One reason the effect of the solar wind flow on reconnection is interesting is that there 
have been conflicting resnlts on whether the cusp reconnection site is stationary or is 
convected tailward. Cowley and Owen [1989] and Gosling et al. [1991] suggested that 
in order to see sunward flow from a reconnection event poleward of the cusp with a 
stationary X-line (where the reconnecting magnetic held goes to zero), the magnetosheath 
how speed should be sub-Alfvenic; if the how is super-Alfvenic, the X-line would have to 
convect tailward. If the magnetosheath how is more than double the magnetosheath 
Alfven speed, reconnection could not occur. Tailward convection of an X-line has been 
seen in global magnetospheric simulations using magnetohydrodynamic simulations of 
Earth [Berchem et al., 1995] and hybrid (kinetic ions with huid electrons) simulations of 
Mercury [Omidi et al, 2006]. In the latter, once the X-line convected far enough tailward, 
a new X-line formed. This behavior was identihed in Cluster observations [Hasegawa 
et al, 2008]. In other studies, the stability of auroral signals associated with high-latitude 
reconnection suggest reconnection sites remain stationary [Fuselier et al, 2000a; Frey 
et al, 2003], although there are uncertainties about whether a lack of change in auroral 
signatures necessarily precludes repeated X-line generation and tailward convection. It 
was also suggested that reconnection with magnetosheath flow speed exceeding twice the 
magnetosheath Alfven speed can occur because suppressing reconnection would introduce 
a pileup of magnetic flux that creates plasma depletion layers, which increase the local 
Alfven speed [Fuselier et al, 2000b]. 
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Another reason this topic is interesting is that a flow shear, such as that caused by 
the solar wind, slows down the reconnection process and can even stop it. It has been 
shown analytically and numerically using the magnetohydrodynamic (MHD) model for 
symmetric systems (with equal and opposite reconnecting magnetic fields and the same 
mass density on either side of the reconnection site) that a super-Alfvenic flow shear 
completely suppresses reconnection, while reconnection still occurs for sub-Alfvenic flow 
shear [Mitchell Jr. and Kan, 1978; Chen and Morrison, 1990; La Belle-Hamer et al., 1994], 
When reconnection occurs with a sub-Alfvenic flow shear present, there is a decrease in 
the reconnection rate [Chen et al., 1997; Li and Ma, 2010; Faganello et al., 2010; Cassak 
and Otto, 2011; Voslion et al, 2011; Zhang et al, 2011; Wu et al, 2013] and outflow 
speed [Cassak, 2011]. The situation is more complicated in more realistic models than 
MHD; there are regimes in the Hall-MHD model in which both tearing (the linear form 
of reconnection) and Kelvin-Helmholtz can be simultaneously linearly unstable [Chacon 
et al, 2003]. 

The suppression of reconnection by a flow shear is potentially of broader importance. It 
was suggested that suppression of reconnection by flow shear limits the length of the X- 
line {i.e., the separator) at the dayside [Borovsky et al, 2013; Borovsky, 2013]. However, 
Komar et al. [2015] pointed out reconnection suppression via flow shear is not expected 
to play a role for southward IMF orientations because the bulk flow is oriented out of the 
reconnection plane (along the flanks) rather than parallel to the reconnecting magnetic 
held (toward the poles), although it does locally decrease the reconnection rate when the 
IMF is directed northward. 
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The effect of flow shear on reconnection is also expected to be relevant at other planets. 
Earth’s solar wind-magnetospheric interaction is qualitatively different from Jupiter’s and 
Saturn’s, where the planet’s relatively rapid rotation contributes to the global convection 
pattern [Vasyliunas, 1983]. Studies have investigated the extent that reconnection occurs 
at the outer planets’ magnetopauses and whether flow shear plays a role in preventing it 
[Masters et al, 2012; Desroche et al, 2012, 2013; Masters, 2014; Fuselier et al, 2014b]. 
The solar wind-magnetospheric interaction at Mercury is similar to Earth’s but on a much 
more rapid time scale [Slavin et al., 2009], so flow shear may also affect reconnection at 
Mercury. 

For applications at Earth’s magnetosphere, it is important to note that the magneto- 
spheric magnetic held is typically a few times stronger than the magnetosheath magnetic 
held and the solar wind plasma in the magnetosheath typically has a much higher density 
than that of the magnetosphere [Phan and Paschmann, 1996; Ku and Sibeck, 1997], i.e., 
the reconnection is asymmetric. Consequently, it is crucial to extend studies of how shear 
to asymmetric systems. There has been much work of late on asymmetric reconnection in 
the absence of how shear; we summarize only those results most germane to the present 
study. The rate of two-dimensional (2D) asymmetric reconnection (with anti-parallel 
magnetic helds) has been studied [Borovsky and Hesse, 2007; Cassak and Shay, 2007]; it 
will be summarized in Sec. 2. In addition, it was found that the X-line and the stagnation 
point, where the inhowing plasma bulk how speed goes to zero, are generally not in the 
same location for asymmetric reconnection [Priest et al., 2000; Siscoe et al, 2002; Dorelli 
et al, 2004; Cassak and Shay, 2007]. 
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We know of only a few nnmerical stndies of the impact of flow shear on asymmetric 
magnetic reconnection. La Belle-Hamer et al. [1995] stndied reconnection with an asym¬ 
metric density and a flow shear in flnid simnlations; they snggested that reconnection is 
snppressed if the flow shear exceeds the Alfven speed on either side of the layer. Tanaka 
et al. [2010] used kinetic particle-in-cell simulations to study reconnection with a density 
asymmetry, a flow shear, and a guide held; they observed that the X-line convects in the 
outhow direction with contributions from both the how shear and the diamagnetic drift 
[Swisdak et al, 2003]. 

In the present study, we use theoretical and numerical techniques to study asymmetric 
reconnection with arbitrary upstream parallel how speeds. We predict the bulk convection 
speed of an isolated X-line using a simple huid analysis. An interesting conclusion is that 
the asymmetries introduce qualitative diherences compared to symmetric reconnection. 
In particular, if the upstream how is equal and opposite on the two sides, the X-line is 
stationary for symmetric reconnection but convects for asymmetric reconnection. We also 
predict the reconnection rate of asymmetric reconnection with upstream how, including 
a condition for the critical upstream how speed required to suppress reconnection. The 
critical how shear (half the diherence of the hows on either side) for symmetric reconnec¬ 
tion is the Alfven speed, so one might expect the cutoh for asymmetric reconnection to 
be the asymmetric generalization of the Alfven speed. However, we show that the cutoh 
exceeds the asymmetric Alfven speed, and the cutoh speed becomes much larger than 
the asymmetric Alfven speed when the asymmetry is large, such as the typical conditions 
at Earth’s magnetopause. Consequently, an isolated X-line at Earth’s polar cusps would 
almost never be suppressed by how shear. We use 2D two-huid simulations to conhrm the 
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predictions. (We point out that particle-in-cell simulations are ostensibly a better tool 
to test the theory because of their more accurate description of kinetic scale physics and 
plasma mixing in collisionless plasmas, but fluid simulations are employed here as a first 
step because they accurately portray the large scale physics while being less noisy and 
therefore easier to compare to the new theory presented here. We employ this approach 
because it has proven prudent for asymmetric reconnection without flow shear [Cassak 
and Shay, 2007, 2008, 2009; Malakit et al., 2010]. Important features which are present 
in particle-in-cell simulations but not two-fluid simulations are discussed in detail in the 
conclusion section.) We then show that the predictions are consistent with recent obser¬ 
vations of a cusp reconnection event for which the magnetosheath flow exceeded twice the 
magnetosheath Alfven speed. The present results are in stark contrast to the previous 
understanding of the effect of flow shear on reconnection based on the Cowley and Owen 
[1989] results, as X-lines convect even for sub-Alfvenic flow and reconnection occurs for 
magnetosheath flow speeds much greater than the magnetosheath Alfven speed. 

The layout of this paper is as follows: Section 2 has the derivation of an expression for the 
convection speed of the X-line and a prediction for the reconnection rate for asymmetric 
reconnection with a flow shear. Section 3 reviews the numerical techniques and parameters 
for the simulations. Section 4 presents the results of our simulations. Section 5 discusses 
implications for observations and applications to planetary magnetospheres. Section 6 
summarizes the results and discusses limitations of the study. 

2. Theory 

We begin by defining system variables for asymmetric reconnection with upstream flow 
parallel or anti-parallel to the reconnecting magnetic helds. The upstream magnetic helds 
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above and below the dissipation region are Bi and B 2 , which are assumed anti-parallel 
[i.e., there is no guide held). The magnetic helds are in opposite directions, so we will use 
Bi and B 2 as the magnitudes of the helds. (Their direction does not impact the present 
analysis as long as they are oppositely directed.) The upstream mass densities are pi and 
p 2 , and the density in the downstream region is pout- 

The upstream how speeds vl,i and vip are dehned in a stationary frame (with the planet 
in question at rest or, in the case of a simulation, the rest frame of the simulation), where 
L refers to the reconnecting magnetic held direction (as in boundary normal coordinates). 
Each speed is dehned as positive if to the right and negative if to the left. The convection 
speed of the X-line is dehned as Udrift- The inhow speeds (normal to the dissipation region) 
are Vinp and fin, 2 , and Ld and 5 are the half-length and half-width of the dissipation region, 
respectively. 

It is convenient to analyze this system in the rest frame of the X-line, so we transform 
into the reference frame moving at a speed Udrift relative to the stationary reference frame. 
The dissipation region in this reference frame is sketched in Fig. 1. In this frame, the 
upstream parallel how speeds are given hj vip— Udrift and vip — i^drift- The outhow speed 
is Uout; which is expected to be the same in both outhow directions in the rest frame of the 
X-line. We dehne x to be the direction of the outhow (L in boundary normal coordinates). 
A number of potentially important ehects are ignored to keep the analysis tractable; these 
are discussed in Sec. 6. 

2.1. Prediction of the X-line Convection Speed 

We hrst consider the convection of an isolated X-line for 2D anti-parallel asymmetric 
reconnection with arbitrary upstream parallel how speeds. The physical reason that the 


DRAFT 


September 7, 2015, 12:29am 


DRAFT 


DOSS ET AL.: ASYMMETRIC RECONNECTION WITH FLOW SHEAR 


X - 11 


X-line can convect in the outflow direction is that the upstream plasmas carry momentum 
in the outflow direction; by conservation of momentum, if the upstream plasma has a non¬ 
zero net momentum, then the dissipation region will too. This is true both for symmetric 
and asymmetric reconnection. For symmetric reconnection, the X-line is expected to 
convect at the average of the upstream flows. For asymmetric reconnection, the inflow 
speed is different for the two upstream sides when the system has asymmetric magnetic 
held strengths [Cassak and Shay, 2007], so the two upstream plasmas do not contribute 
momentum in the outflow direction equally. Consequently, as we will show, the X-line 
convects in the outflow direction for asymmetric reconnection even for equal and opposite 
upstream flow speeds. An alternate, but equivalent, interpretation is that the X-line and 
stagnation point are separated during asymmetric reconnection [Cassak and Shay, 2007], 
so the side with the plasma that crosses the X-line imposes a flow in the outflow direction 
at the X-line, causing it to convect. 

To estimate the convection speed, we use a fluid description of the plasma to And the 
bulk effects of the flow. The governing equation is the momentum equation of MHD, 
which in conservative form is 



( 1 ) 


where p is the plasma density, v is the bulk velocity, P is the gas pressure, B is the 
magnetic held, and I is the unit tensor. Take the volume integral of Eq. (1) over the 
entire dissipation region. For steady-state reconnection, the time derivative term vanishes. 
From the divergence theorem, the remaining term becomes 
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where dS is a differential area element directed normal to the boundary of the dissipation 
region. For both the upstream and downstream boundaries, the magnetic field is approx¬ 
imately parallel to the boundaries, so the BB term does not contribute. Total pressure is 
balanced across the dissipation region, so the pressure term has no net contribution. The 
X component of the surviving term of Eq. (2) is 

y dS ■ (pvv^) = 0. (3) 

Performing a scaling analysis on this gives 




( 4 ) 


where the two terms are the contributions from the two upstream sides and the x directed 
momentum flux on the downstream edges cancel. Conservation of magnetic flux implies 
VinpBi ~ nm, 2 -B 2 [Cassttk and Shay, 2007], so solving for Udrift gives 


"^drift 


PiB2Vl,i + P2BiVl,2 
P 1 B 2 + P 2 BI 


( 5 ) 


This gives the convection speed of the X-line in the outflow direction for arbitrary up¬ 
stream densities, reconnecting magnetic field strengths, and parallel flow speeds. As 
discussed earlier, Udrift is non-zero for symmetric anti-parallel flow {vl^i = —ung) when 
the magnetic fields and/or densities are asymmetric. This shows that the X-line convects 
regardless of upstream flow speeds, not only for super Alfvenic flow as suggested in Cowley 
and Owen [1989]. 

We consider particular limits of this expression. For any upstream conditions, if = 
vl ,2 = 't'shear, then Udrift = 't'shear, as expected. For the case of symmetric reconnection with 
Bi = B 2 and pi = p 2 , Eq. (5) reduces to the expected result of 


Wrift 


Vl,1 + Vl,2 


( 6 ) 
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the average of the two upstream flow velocities (which is zero if the flows are equal and 
opposite). If the densities are symmetric but the magnetic helds are not, then 


'Adrift 


B2VL,1 + BiVl,2 

B2 + Bi 


and if the magnetic helds are symmetric but densities are not, we hnd 


■ydrift 


PiVL,1 + P2VL,2 
Pl+ P2 


These predictions are testable with simulations. 


(7) 


( 8 ) 


2.2. Prediction of the Reconnection Rate 

We turn to the prediction of the reconnection rate of 2D asymmetric anti-parallel re¬ 
connection with arbitrary upstream parallel how speeds. The reconnection rate T'shear.sym 
for symmetric reconnection with an equal and opposite upstream how was recently found 
to scale as [Cassak and Otto, 2011] 


E, 


shear,sym 


En ( 1 


•^shear 


(9) 


where Eq is the reconnection rate in the absence of upstream how (~ 0.1 for collisionless 
reconnection), fshear is the upstream speed of the symmetric plasma bulk how, and ca is 
the outhow speed given by the Alfven speed based on the reconnecting magnetic held. 
Physically, the form of this correction is motivated in a way similar to the suppression 
of reconnection by diamagnetic ehects. Swisdak et al. [2010] argued that the magnetic 
tension force on the plasma due to a newly reconnected magnetic held line has to overcome 
the momentum of the moving plasma in the dissipation region. 

For asymmetric reconnection, we propose that the ehect is very similar. One main dif¬ 
ference is that the outhow speed during 2D anti-parallel asymmetric reconnection becomes 
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[Cassak and Shay, 2007; Swisdak and Drake, 2007] 

2 B1B2 Bi + B2 

CA,a.ym ~ 47r ^ 1^2 + ^2^1 ’ 

A second difference is that the X-line and stagnation point are not at the center of the 
dissipation region [Cassak and Shay, 2007], so the upstream plasmas on the two sides 
contribute momentum in proportion to the relative location of the stagnation point. This 
is illustrated in Fig. 2, showing a recently reconnected magnetic field line in blue with 
the stagnation point denoted by the “S”. The distance from the upstream edges to the 
stagnation point above and below are and 632 , and the left side shows how the two 
upstream flows impact the straightening magnetic field line. Thus, not only is the up¬ 
stream momentum allowed to be different on the two sides because the density and flow 
speed can differ, but they contribute in different proportions to the dissipation region. 

We follow the result of Eq. (A3) from Swisdak et al. [2010] and propose the outflow 
speed decreases as 



^2 

^A,asym 


ox (l’L,! L’drift) „ ^ (^T,2 L’drift) 


52 , 


25 


25 


( 11 ) 


In writing this, we effectively subtract off the energy due to the flow, but here we do so 
in proportion to each side’s contribution opposing the release of magnetic energy. Using 
5si/25 = P 1 B 2 /{P 1 B 2 + P 2 B 1 ) and 5s2/2(5 = p 2 -Bi/(pi-B 2 + P 2 -B 1 ) [Cassak and Shay, 2007] 
and eliminating Udrift using Eq. (5), some algebra reveals 


~ 4,asym “ {^L,! “ VLpf 


P1B2P2B1 

{P1B2 + p 2 Biy’ 


( 12 ) 


Note, the only upstream flow dependence is on the difference in upstream flows vl^i — ul, 2 - 
Consequently, we suggest the reconnection rate Eshear.asym for asymmetric reconnection 
with arbitrary upstream parallel flow speeds, extending Eq. (9) to asymmetric systems. 
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scales as 


-^shear,asym -^0,asym I 1 


■^shear 


4pii?2P2-Bl 


where 


^0,asym 


Ci,asym (pA + ’ 

B\B2 Cyi^^asym 2(5 


(13) 


(14) 


Bi + i ?2 c Ld 

is the asymmetric reconnection rate in the absence of upstream flow [Cassak and Shay, 
2007] and 


^shear 


VLA - Vt2 


(15) 


is half the difference of the upstream flow speeds (their average if oppositely directed). 
Thus, the reconnection rate decreases with increasing flow shear. 

An important consequence of Eq. (13) is the prediction of a critical flow speed Ushear.crit 
above which reconnection does not occur (corresponding to Eshear.asym = 0). The predic¬ 
tion is that reconnection shuts off above flow shear speeds of 


P 1 B 2 + P 2 B 1 

^shear.crit asym":r 7 ~ rTTTTo * 

2{piBip2BiYI‘ 


(16) 


In the symmetric reconnection limit, this expression reduces to Ushear.crit ~ ca, the known 
result, so it generalizes the known result to asymmetric systems. It is important to note 
that the fraction multiplying CA,asym is always greater than or equal to one. This im¬ 
plies that the critical flow shear required to suppress asymmetric reconnection exceeds 
the asymmetric Alfven speed. For larger asymmetries, the critical flow shear becomes 
larger. Consequently, while super-Alfvenic flow is sufficient to suppress symmetric re¬ 
connection, reconnection can proceed in asymmetric reconnection for super-Alfvenic flow, 
which differs from the Cowley and Owen [1989] prediction. If the asymmetry is large (as 


DRAFT 


September 7, 2015, 12:29am 


DRAFT 






X - 16 


DOSS ET AL.: ASYMMETRIC RECONNECTION WITH FLOW SHEAR 


is typical at Earth’s magnetopause), then the suppression condition is much larger than 
the asymmetric Alfven speed. This is discussed further in Sec. 5. 

3. Simulation Setup 

To test the predictions, we perform 2D simulations using the massively parallel two- 
fluid code F3D [Shay et al, 2004]. The code solves for the density, ion velocity, magnetic 
held, and ion pressure (assumed to be an adiabatic ideal gas with ratio of specihc heats 
7 = 5/3). Electrons are assumed cold for simplicity. The numerical algorithm is the 
trapezoidal leapfrog in time and fourth order hnite difference in space. 

The code evolves variables normalized to Bq, the reference magnetic held strength, and 
pQ, the reference mass density, where these quantities typically are the initial values on 
one upstream side of the simulation. Other variable’s normalizations are derived from 
these values: velocities are normalized to the Alfven speed cao = Bo/{4:7rpoY^^, lengths 
are normalized to the ion inertial length dio = {rriiC^/, times are normalized to 
the inverse ion cyclotron frequency 0“^ = {eBQ/mic)~^, electric helds are normalized to 
caoBq/c, and pressures are normalized to Bq/Att, where ruj and e are the ion mass and 
charge and uq = po/rrii is the plasma number density. The x, y and 2 ; directions are aligned 
with the initial direction of the magnetic held (L in boundary normal coordinates), the 
inhow (A^), and the out-of-plane current (M). 

Periodic boundary conditions are used in each direction. The size of the computational 
domain is L^x Ly = 204.8 x 102.4 dio, and the grid scale is 0.05 dio in each direction. The 
electron mass is me = mi/2b, and it is not expected that the large-scale behavior seen in 
our simulations is sensitive to this value. 
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The initial conditions for the simulation have a magnetic held prohle given by an asym¬ 
metric double Harris sheet, 



Ly/4 < \y\ < Lyl2 
0 < |i/| < Ly/4 


(17) 


with an initial current sheet width of = dio. We do not use an initial out-of-plane 
(guide) magnetic held. The initial density prohle is 



(18) 


with asymptotic values poi in the central portion of the domain and po 2 at the top and 
bottom of the domain. The initial pressure prohle is chosen to balance the pressure 
identically, with a minimum value of where = Hiax(i?oi, H 02 ) is the 

stronger of the two magnetic held strengths and (dram is the minimum plasma beta on 
either side of the plasma. We choose d min = 2. 

The bulk how is initialized with a prohle of 



(19) 


We use the same tco as for the magnetic held prohle (although this is not a requirement 
[Li and Ma, 2010]). The upstream how speeds are vi and ^ 2 , which for the simulations 
presented here, are always equal and opposite: Vi = —V 2 - We dehne this common speed 


as Ushean which poiuts to the left at the top and bottom of the computational domain and 


to the right in the central part of the domain. 

The simulations do not employ any explicit resistivity or viscosity. However, there is 
a fourth order dihusion used in each evolution equation to damp noise at the grid scale. 


with coefhcient 1 x 10 ^ for asymmetric magnetic held simulations and 1.25 x 10 ® for 


asymmetric density simulations. The magnetic held is perturbed initially with random 
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noise of amplitude 0.00005 Bq to break symmetry, which allows any secondary magnetic 
islands that arise to be ejected. A coherent magnetic perturbation of the form (5B = 
— (0.012i?of^y/27r)z X V[sm.{2T{x/L^) sit ?{21111 jLy)] is used to initiate reconnection in a 
controllable manner. 

A series of numerical simulations is performed with initial magnetic asymmetries of 
3 or 2 with uniform density and with plasma density asymmetries of 3 with symmetric 
magnetic held. The initial how speed Ughear in the upstream regions is varied between 0 
and 2.4 c^o- Table 1 shows a list of the various simulations from which data are collected. 

4. Results 

We begin with a case study showing that the X-line in asymmetric reconnection convects 
in the outhow direction for equal and opposite upstream how speeds. The simulation has 
asymmetric magnetic helds Bi = 3 Bq and B 2 = Bq and symmetric density pi = P 2 = Po 
with an upstream how speed of Ughear = 1-2 cao- Figure 3 shows the out-of-plane current 
density as the background, with magnetic held lines overplotted in blue for a portion 
of the computational domain. The top, middle, and lower panels are at times t = 90,105, 
and 120 respectively. The X-line, located near the peak in the current density, clearly 
convects to the right as shown by the white arrow. It travels approximately 18 dio in 30 
giving a convection speed close to 0.6 cao, which is the predicted speed from Eq. (5). 

To extract more quantitative information from the simulations, the systems are evolved 
until a steady-state is reached with the X-line convection speeds and reconnection rates 
being relatively constant. To apply a common criterion for the steady-state across sim¬ 
ulations (which have diherent reconnection rates), we put a condition on the size of the 
primary magnetic islands; all simulations presented here are relatively steady when the 
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island half widths are between 7 and 11 dio. X-line convection speeds and reconnection 
rates are measured during these times. 

The X-line and 0-line are found in the standard way as the saddle point and local 
extremum, respectively, near the current sheets of the flux function ip given by B = z x Vip. 
The X-line convection speed is obtained from the time average of the time derivative of 
the position of the X-line during the steady-state interval. The reconnection rate is the 
average of the time rate of change of the difference in magnetic flux between the X-line 
and the 0-line. 

4.1. The X-line Convection Speed 

For the X-line convection speed, we hrst use simulations with asymmetric magnetic 
helds, symmetric plasma densities, and symmetric upstream flow. The x coordinate of 
the X-line in the simulation reference frame is plotted as a function of time in Fig. 4. 
The (blue) solid line is from a simulation with i?i = 3 i?o, B 2 = Bq and the (red) dashed 
line is for Bi = 2 Bq, B 2 = Bq, both with Ushear = 1-2 caq. The convection speed is higher 
with the stronger magnetic held, as predicted, with measured values of Udrift = 0.63 cao 
for Bi = 3 Bo and 0.37 cao for Bi = 2 Bq, compared to predicted values from Eq. (7) of 
0.60 and 0.40 cao- 

To test the upstream how speed dependence. Fig. 5 shows the average drift speed Udrift 
as a function of upstream how speed Ushear for simulations with Bi = 3 Bq and B 2 = Bq. 
Here, and throughout, the (blue) triangles are for the top current sheet and the (red) 
squares are for the bottom current sheet. The prediction from Eq. (7) is plotted as the 
dashed line, and the results clearly agree well with the prediction. 
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Simulations with symmetric magnetic fields but with asymmetric plasma densities were 
also carried out. However, we do not expect the numerical results to be reliable because 
of a known problem with the fluid approach in systems with asymmetric density. In 
particular, fluid simulations require conduction to allow for plasma mixing [Cassak and 
Shay, 2009]; newly reconnected field lines have different densities and temperatures, and 
in the absence of conduction the strong parallel temperature gradient persists instead 
of the plasmas mixing. These simulations do not contain conduction and it is unlikely 
the standard fluid closure reproduces the more realistic kinetic mixing in a collisionless 
dissipation region. Particle-in-cell simulations will be required to assess this prediction 
for asymmetric densities. 

4.2. The Reconnection Rate 

For the reconnection rate in systems with asymmetries and upstream flow, we plot the 
reconnection rate E as a. function of flow speed Ushear in Fig. 6(a) for simulations with 
fixed Bi = 3 Bq and B 2 = Bq for upstream flow speeds from 0 to 1.6 cao- The prediction 
from Eq. (13) is plotted as the dashed line, where Eo^asym is chosen as the measured 
value of the reconnection rate in the simulation with no upstream flow. The agreement 
is very good. The result for Ushear = 1.2 cao falls a little below the curve; this simulation 
had a significant secondary island unlike the others which motivates its slightly worse 
performance. 

A similar study is carried out with simulations of asymmetric plasma densities but with 
symmetric magnetic field. While the convection speed is not expected to be correct in 
these simulations, it was argued that the reconnection rate is reliable [Cassak and Shay, 
2009; Birn et ai, 2010]. (This is because redoing the calculation of the outflow speed 
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and reconnection rate from Cassak and Shay [2007] for a system where the downstream 
plasma does not mix leads to the same expressions as for the system where mixing occurs. 
The differences arise only in the sub-structure of the dissipation region.) The results for 
a simulation study with pi = po and p2 = 3 Po with upstream flow speeds varied from 0 
to 0.6 cao are shown in Fig. 6(b). The prediction from Eq. (13) is plotted as the dashed 
line. For asymmetric densities, the agreement with the prediction is again very good. 

4.3. The Cessation Condition 

We now test the condition for the cessation of reconnection, given by Eq. (16). For 
the series of simulations with Bi = 3 Bq,B2 = Bq, pi = P2 = Po, we have c^,asym = 
\/3 Cao from Eq. (10), so Ushear.crit = 2 Cao- In the simulations, we see that reconnection 
occurs, with the predicted reconnection rate and X-line drift speed, for Ughear = 1-6 cao- 
In contrast, for Ushear = 2.4 cao, there is a phase change and reconnection is not the 
dominant effect. This is shown in Fig. 7, which gives the out-of-plane current density 
Jz for this simulation. There is clear evidence of the early phases of a Kelvin-Helmholtz 
instability. We note that reconnection is known to occur as a secondary instability of 
Kelvin-Helmholtz in a plasma, and that does occur here, but this is fundamentally different 
than the reconnection-dominated current sheets displayed in Fig. 3. We note that for 
Vsheur = 2 cao, fight at the predicted cutoff, there is a hybrid of reconnection and Kelvin- 
Helmholtz rolls forming. Thus, the results are consistent with the prediction. As the 
qualitative behavior near the cutoff becomes more challenging to assess, it is prohibitive 
to pin down the transition flow speed to higher precision than done here. 

The asymmetric density simulations with Bi = B 2 = Bq, pi = pQ, and p 2 = 3 po give 
another opportunity to test the prediction. Here, c^,asym = \/l/2 cao, so nshear,crit = 
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^2/3 caq — 0.82 caq- In the simulations, reconnection does occur for Ushear = 0.6 cao but 
is not appreciable for Ushear = 0.8 cao, which is consistent with the predictions. 

In summary, the predictions for Udrift in Eq. (5), Eshear,asym in Eq. (13), and Ushear,crit in 
Eq. (16) agree with the simulation results. 

5. Applications to Planetary Magnetospheres 

For a sample application, we consider reconnection at the cusps of Earth’s magneto¬ 
sphere. We emphasize that the predictions in Sec. 2 assume an isolated X-line; we hrst 
discuss how the present results suggest isolated X-lines would act in the magnetosphere. 
We call side 1 “ms” for magnetosphere and side 2 “sh” for magnetosheath. For typical 
conditions at the magnetopause, psh S> Pms and > Bgh- Also, the upstream flow 
speed in the magnetosphere VL,ms is negligible, while the magnetosheath flow VL,sh is due 
to the solar wind flowing around the magnetopause. In this limit, the X-line convection 
speed from Eq. (5) becomes 

Wrift — '^L,sh- (2b) 

Thus, an isolated X-line would convect tailward at essentially the same speed as the flow 
in the magnetosheath, not the average flow speed as one would expect for symmetric 
reconnection. The physical cause for this, as sketched in Fig. 2, is that the stagnation 
point is far to the magnetospheric side of the dissipation region, so most of the dissipation 
region is populated by magnetosheath plasma and most of its momentum is contained 
in the magnetosheath plasma. Therefore, it flows near the speed of the magnetosheath 
plasma. 

This is borne out in recent Cluster observations of a tailward convecting X-line [ Wilder 
et al., 2014] which we use as a case study. In this event at the southern hemisphere’s cusp 
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on December 27th, 2005, Cluster’s Cl and C3 spacecraft witnessed reconnection signatures 
on a crossing from the magnetosheath to the magnetosphere. The magnetosheath had 


Bgh — 10 — 15 uT and a number density Ugh ~ 60 — 70 cm while in the magnetosphere 


Bms — 60 nT and Ums = 0.5 cm“^. From their separation and the time delay between 
observed jet reversals, the convection speed of the X-line was estimated to be Udrift = 105 
km/s, while the L component of the magnetosheath flow was estimated at VL,sh = 106 
km/s. This is consistent with Eq. (20). 

In the past, assessing whether reconnection could occur was done by comparing the 
magnetosheath flow speed to the magnetosheath Alfven speed [Cowley and Owen, 1989]. 
For this event, the magnetosheath Alfven speed is 28 km/s, while the magnetosheath flow 
speed is 105 km/s. This exceeds twice the magnetosheath Alfven speed so previous models 
would suggest reconnection should not occur, but reconnection is observed to happen. The 
asymmetric Alfven speed for these parameters is 74.5 km/s, so Ushear = i^i,,sh/2 — 53 km/s 
is sub-Alfvenic, and the prediction here is that reconnection should occur. 

For the reconnection rate, in the limit appropriate for typical conditions at the magne¬ 
topause, Eq. (13) becomes 



( 21 ) 


Since psh S> Pms-, this implies the surprising result that the reconnection rate is not 
changed much by the flow shear for typical magnetospheric parameters. This is a major 
departure from the symmetric case where the reconnection rate falls strongly as a function 
of flow shear speed. Physically, the cause for this is related to the stagnation point being 
very close to the magnetospheric side, as before. In the reference frame of the X-line, the 
magnetosheath is essentially still, while the magnetosphere has a signihcant flow. However, 
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since the magnetospheric density is so low, it has very little effect on the reconnection 
process. 

This is seen fnrther by investigating the condition for snppressing reconnection. From 
Eq. (21) and nsing Wshear — 'yL,s/i/2 from Eq. (15) since the magnetospheric plasma is 
essentially stationary, the condition that reconnection is suppressed (Eghear,asym < 0) is 

^ f PshBms\ /nn\ 

'^LjSh ^ ^AjSsym I ^ 1 • (22) 

\ Pms-Dsh J 

For the example of the case study event from Cluster, this implies it would take a mag¬ 
netosheath flow speed of 22 times bigger than CA,asym — [Bsh{Bms + Bsh)/^T^Psh\^^‘^ for 

flow shear to suppress reconnection, far greater than in symmetric reconnection and the 
Cowley and Owen [1989] model where Alfvenic flow suppresses reconnection. (We point 
out further that particle-in-cell simulations and observations typically hnd outflow speeds 
half of that predicted in symmetric or asymmetric reconnection theories, so the critical 
speed may be a factor of two smaller than this prediction. However, this does not change 
the result that the critical speed is much higher than the asymmetric Alfven speed and the 
magnetosheath Alfven speed.) Thus, for an isolated X-line with typical magnetospheric 
parameters, the flow shear would rarely prevent reconnection from occurring. 

While X-lines in the magnetosphere can convect as discussed earlier, it was suggested 
from observations [Fuselier et ai, 2000a; Frey et al, 2003] that cusp reconnection can be 
steady for hours with a stationary X-line, though it may be difficult to be sure solely from 
auroral signatures. Further, the distance to the X-line has been inferred from observations 
[Petrinec et al, 2003; Trattner et al, 2007a, b], and strong convection is not seen, even 
though Eq. (20) would predict a large enough convection speed to be observable. This 
suggests that reconnection in the magnetosphere is not always consistent with isolated 
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X-lines, which was a key assumption of the theory presented here. In particular, the 
theory does not take any external effect other than flow into account. A single dayside 
X-line is not isolated; instead, for a single X-line during steady times to be stationary, it 
must balance all magnetic forces, gas pressure forces, and dynamic pressure forces from 
bulk flow (of the type discussed here). Thus, by definition, the flow effects have already 
been taken into account to set up a single X-line, and the present model would not apply. 
Physically, the present model does not include magnetic helds in the cusp being line-tied 
to the ionosphere, which is an additional force that can slow or prevent X-line convection 
and would undoubtedly modify the present results. It is possible that the tailward moving 
X-lines seen in observations [Hasegawa et ai, 2008; Wilder et ai, 2014] are related to a 
secondary island moving tailward; a secondary X-line is not in the same equilibrium as 
a single X-line, so it can convect tailward. Future work is necessary to determine the 
conditions under which some X-lines convect and others do not, the speed at which they 
convect, and their signatures in the cusp. 

A second point worth making is that the results in this section would change in the 
presence of a high-density plasmaspheric drainage plume [Borovsky and Steinberg, 2006; 
Borovsky and Denton, 2006] on the magnetospheric side of the dissipation region. This 
increases the local density, which decreases the outflow speed [Walsh et ai, 2014a, b]. 
Plumes are most commonly observed centered around 13.6 h magnetic local time [Walsh 
et ai, 2013]. If the magnetospheric density is comparable or larger than the magne¬ 
tosheath, the X-line would not convect with the magnetosheath flow, and the full expres¬ 
sions from Sec. 2 would be needed. 
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We note in passing that the signihcant departures between symmetric and asymmetric 
reconnection with a flow shear suggest that the effect of flow shear at Jupiter, Saturn, 
and Uranus should be revisited. The previous work on this subject employed the predic¬ 
tions for symmetric reconnection, so the results are likely to change with the asymmetric 
reconnection results presented here. 

It is worth putting the present results in context of a leading model for what allows 
cusp reconnection to occur despite nominally super-Alfvenic flow. It was argued that 
the draped magnetic held during northward IMF compresses, and therefore increases, the 
magnetopause magnetic held, which itself causes the mass density to decrease, forming a 
density depletion layer at the magnetopause [Petrinec et ai, 2003]. Both the increased 
held and decreased density increase the local Alfven speed, potentially making the flow 
sub-Alfvenic. The present theory should not be interpreted as implying that a density 
depletion layer does not occur or that the layer would not change the local Alfven speed. 
Using the Wilder et al. [2014] parameters, Petrinec et al. [2003] would suggest (assuming 
that the Alfven Mach number is based on magnetosheath parameters) that a density 
depletion by a factor of about 2.5 would be sufficient to make the flow sub-Alfvenic (see 
their Fig. 7). When based on the asymmetric Alfven speed, the necessary depletion 
factor reduces to approximately 1.3. Clearly, more research is needed to identify the level 
of plasma depletion, its effect on cusp reconnection, and the extent to which the present 
results apply to the dayside magnetopause. 

6. Discussion 

In this study, we use a scaling argument based on conservation of momentum to And the 
convection (drift) speed of the X-line in anti-parallel 2D asymmetric magnetic reconnec- 
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tion with arbitrary upstream parallel flow speeds [Eq. (5)]. We also present a prediction 
for the reconnection rate for arbitrary upstream flow speeds [Eq. (13)]. The predictions 
are confirmed using 2D two-fluid numerical simulations. For asymmetric magnetic fields, 
the results agree well for both predictions. For asymmetric densities, simulations agree 
with the reconnection rate prediction, but the drift speed cannot be assessed because 
the fluid model does not correctly model plasma mixing. In particular, the results show 
that the X-line convects even for sub-Alfvenic flow, which contrasts the leading model by 
Cowley and Owen [1989]. 

The reconnection rate prediction gives a threshold flow shear above which reconnec¬ 
tion does not occur [Eq. (16)], and this prediction is consistent with the simulations. 
This result shows that asymmetric reconnection can persist with flow shear exceeding the 
asymmetric Alfven speed, and can occur even with much larger flow shear speeds if the 
asymmetries are large. We note that the critical speed prediction in Eq. (16) differs from 
the condition given by La Belle-Hamer et al. [1995], who suggested the critical speed is the 
larger of the two Alfven speeds on either upstream side of the dissipation region. It also 
differs from the prediction that reconnection does not occur if the difference between the 
magnetosheath and magnetospheric flow exceeds twice the magnetosheath Alfven speed 
[Cowley and Owen, 1989]. To see this, consider the simulation with Bi = 3,i?2 = 1, a 
uniform density of pi = p 2 = 1, and a flow of Ushear = 1-6. Here, the magnetosheath nom¬ 
inally corresponds to the “2” side (with the weaker magnetic held). In the rest frame of 
the magnetosphere (the “1” side), the magnetosheath flow speed is 3.2. This exceeds dou¬ 
ble the magnetosheath Alfven speed of 1. The fact that reconnection is observed in this 
simulation is evidence against the Cowley and Owen [1989] model. The fact that we do 
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not see reconnection in the fshear = 2.4 simnlation is evidence against the La Belle-Hamer 
et al. [1995] model, which is below the magnetospheric Alfven speed of 3. 

The results have potentially important implications for reconnection at Earth’s magne¬ 
topause. For isolated X-lines, the predictions suggest that the X-line convects essentially 
with the magnetosheath flow as a consequence of the stagnation point being nearly all 
the way to the magnetospheric side of the dissipation region. The reconnection rate is 
affected by only a small amount, and it would take magnetosheath flow an order of mag¬ 
nitude faster than the asymmetric Alfven speed to suppress reconnection. This is a major 
departure from the current understanding based on symmetric reconnection which claims 
magnetosheath flow greater than twice the Alfven speed is sufficient to suppress reconnec¬ 
tion. The results may differ for non-isolated X-lines, such as magnetic helds line-tied to 
the ionosphere during reconnection near the polar cusps. The results could also drastically 
alter previous estimates of how flow shear affects dayside reconnection at Jupiter, Saturn, 
and Uranus, so revisiting these results is important future work. Another potential ap¬ 
plication is determining whether the present results impact predictions for the speed of 
anti-sunward propagation of flux transfer events [Cowley and Owen, 1989; Cooling et al, 
2001 ]. 

The present analysis made a number of simplifying assumptions that should be ad¬ 
dressed in future work. Since the analytical calculation uses the fluid picture, hnite 
Larmor radius effects are ignored, but they may be important in the boundary layers 
[Malakit et al, 2013; Koga et al, 2014], especially for the dayside magnetopause where 
the density asymmetry is typically signihcant. The present analysis and simulations ig¬ 
nore asymmetries in the outflow direction [Murphy et al, 2010; Oka et al, 2011], which 
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may be important for reconnection near the polar cusps. In particular, asymmetries in 
the outflow direction cause the outflow speeds to be different in the two outflow directions 
[Murphy et al, 2010], so future work is necessary to determine the effect of shear flow 
in such systems. The present analysis and simulations also do not include a guide held, 
which may be present [Muzamil et ai, 2014]. The guide held is important because, when 
coupled with a gas pressure gradient across the inhow direction at the dissipation region, 
it can set up diamagnetic ehects which also have been shown to cause the X-line to convect 
in the outhow direction [Swisdak et al, 2003; Beidler and Cassak, 2011], and therefore 
either reinforce or oppose the convection caused by how shear [Tanaka et al, 2010]. The 
analysis and simulations also do not include upstream hows in the out-of-plane direction 
[Wang et al, 2012; Chen et al, 2013; Wang et al, 2014; Tassi et al, 2014]. Also, the simu¬ 
lations used the huid model, which does not self-consistently capture plasma mixing in the 
exhaust for systems with asymmetric density [Cassak and Shay, 2009]; this would be hxed 
using simulations employing the particle-in-cell technique [Roytershteyn and Daughton, 
2008; Tanaka et al, 2010]. Another future study should address asymmetric reconnection 
with a how shear in line-tied systems relevant to polar cusp applications. 
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Figure 1. Sketch of a reconnection region for asymmetric reconnection with flow shear 
in the reference frame the dissipation region. Red dashed lines show components of plasma 


velocities and blue solid lines represent the reconnecting magnetic helds. 
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Figure 2. Sketch of a newly reconnected magnetic field line (blue curve) for a system 
with upstream flow and asymmetries. The box denotes the left half of a dissipation region. 
“X” and “S” mark the X-line and stagnation point, and the dashed line goes through the 
stagnation point. How the bulk flow impacts the dissipation region is sketched on the left 
side. Nominally, the relative positions of the X-line and stagnation point are representative 
of the magnetosphere being at the top and the magnetosheath being at the bottom. 
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Table 1. List of simulations with their initial values of upstream magnetic fields Bi and 
B 2 (in units of i?o), mass densities pi and p 2 (in units of po), and upstream flow speed Ushear 
(in units of cao)- Also included is the predicted convection speed Udrift,pred from Eq. (5), 
measured convection speeds of the top (udrift,T) and bottom (udrift,^) current sheets, the 
predicted reconnection rate Epred from Eq. (13) (in units of caqBq/c), and the measured 
reconnection rates from the top Et and bottom Eb current sheets. The Eo,asym value is 
obtained from the average measured reconnection rate in runs without shear, except for 
the Bi = 2 Bq simulation which uses the prediction from Cassak and Shay [2007] using 
5/L = 0.06. 
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Figure 3. Out-of-plane current density as a function of x and y (in units of 
dio), with magnetic field lines overplotted in blue, for a simulation with Ushear = 1-2 c^o, 
i?i = 3 Bo,B2 = Bq, and pi = p2 = Po- The top, middle, and bottom plots are at time 
t = 90,105, and 120 0“^. Only a portion of the computational domain is plotted. 
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Figure 4. X-line position (in units of dio) as a function of time (in units of 
for simulations with Ughear = 1-2 cao for Bi = 3 Bq,B 2 = Bq (blue solid) and Bi = 
2 Bq, i ?2 = Bq (red dashed). The average drift speeds are Udrift = 0.63 cao for Bi = 3 Bq 
and Udrift = 0.37 c^io for Bi = 2 Bq. 
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Figure 5. X-line convection speed Udrift (in units of cao) for the top (A) and bottom 
(□) current sheets from simulations with magnetic helds Bi = 3 Bq and B 2 = Bq a.s a. 
function of upstream flow speed Ughear (in units of c^o)- The dashed line is the prediction 
from Eq. (7). 
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Figure 6. Reconnection rate E (in units of caqBq/c) for the top (A) and bottom (□) 
current sheets from simulations with various flow speed Ushear (in units of cao) and (a) 
magnetic field asymmetry Bi = 3 Bq and B2 = Bq and (b) density asymmetry pi = po 
and p 2 = 3 po. The dashed lines are the predicted rates from Eq. (13), with Fo^asym as 
the average rate of reconnection in the absence of a sheared flow. 
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Figure 7. Out-of-plane current density Jz, with magnetic held lines overplotted in blue, 
for a simulation with Bi = 3 Bo, B 2 = Bq, pi = p 2 = Po with Ushear = 2.4 cao , which is 
above the predicted cutoff where reconnection is suppressed. 
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